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Abstract 



We present a complete set of the equations and matching condi- 
tions required for the description of physicahy meaningful charged, 
dissipative, spherically symmetric gravitational collapse with shear. 
Dissipation is described with both free-streaming and diffusion ap- 
proximations. The effects of viscosity are also taken into account. 
The roles of different terms in the dynamical equation are analyzed 
in detail. The dynamical equation is coupled to a causal transport 
equation in the context of Israel-Stewart theory. The decrease of the 
inertial mass density of the fluid, by a factor which depends on its 
internal thermodynamic state, is reobtained, with the viscosity terms 
included. In accordance with the equivalence principle, the same de- 
crease factor is obtained for the gravitational force term. The effect 
of the electric charge on the relation between the Weyl tensor and the 
inhomogeneity of energy density is discussed. 
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1 Introduction 



The study of self-gravitating spherically symmmetric charged fluid distribu- 
tions has a long and a venerable history, starting with Rosseland and Ed- 
dington's contributions [HE]- Since then a large number of works have been 
dedicated to making manifest the influence of electric charge on the structure 
and evolution of self-gravitating systems (see [lJ-[42] and references therein). 

Although some of these works refer to static situations ([3]- [5], [9l [191 ISSl 
122 [291 [30], I32]-I3l], [361 [371 [391 m]) there have been important efforts in 
describing dynamical situations too ([6]- [8], [I0]-[T1], [Ml [TTl [211 [22l [261 [iQl 
[12]). Particularly relevant for the present paper are references [HI [22] and 

m- 

A renewed interest in this subject emerges from the appearance of new 
mechanisms allowing for the presence of huge electric charge in self-gravitating 
systems. From simple classical considerations, it can be shown that physical 
objects with large amounts of charge (much larger than 100 Coulomb per 
solar mass) cannot exist [H [21 131] • Furthermore, as shown by Bekenstein 
IH] the electric charge is bounded by the fact that the resulting electric field 
should not exceed the critical field for pair creation, 10^^ V cm~^. However, 
these restrictions have been questioned by several authors [151 [HI [23 [37] . 
Particularly appealing is the possibility of very high electric fields in strange 
stars with quark matter (see [5H1 [39J and references therein). 

All this having been said, it should be clear that the restrictions men- 
tioned above refer to equilibrium (stable) configurations. They do not apply 
to phases of intense dynamical activity with time scales of the order of (or 
even smaller than) the hydrostatic time scale, and for which the quasi-static 
approximation is clearly not reliable (e.g. the collapse of very massive stars 
[43j or the quick collapse phase preceding neutron star formation, see for 
example [44j and references therein). The description of this very dynamic 
regime is the main purpose of this manuscript. 

Besides electric charge, which will be assumed to comove with the fluid, 
we shall also consider dissipative phenomena. It is already an established fact 
that gravitational collapse is a highly dissipative process (see [1S]-[1Z] and 
references therein), so the relevance of dissipation in its study cannot be over- 
emphasized. Dissipation due to the emission of massless particles, photons, 
and/or neutrinos is a characteristic process in the evolution of massive stars. 
In fact, it seems that the only plausible mechanism for carrying away the 
bulk of the binding energy of a star collapsing to a neutron star or black hole 
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is neutrino emission |48j . 

In the diffusion approximation, it is assumed that the energy ffux of 
radiation, hke that of thermal conduction, is proportional to the gradient of 
temperature. This assumption is in general very sensible, since the mean free 
path of particles responsible for the propagation of energy in stellar interiors 
is normally very small compared with the typical length of the object. Thus, 
for a main sequence star such as the sun, the mean free path of photons 
at the centre is of the order of 2 cm. Also, the mean free path of trapped 
neutrinos in compact cores of densities above about 10^^ g cm~'^ becomes 
smaller than the size of the stellar core ^91 [50] . 

Furthermore, the observational data collected from supernova 1987A in- 
dicates that the regime of radiation transport prevailing during the emission 
process is closer to the diffusion approximation than to the free streaming 
limit [51j . 

However, in many other circumstances the mean free path of particles 
transporting energy may be large enough to justify the free streaming ap- 
proximation. Therefore it is advisable to include simultaneously both limit- 
ing cases of radiative transport, diffusion and free streaming, allowing us, by 
taking both in combination, to describe a wide range of situations. 

The effects of dissipation, in both limiting cases of radiative transport, 
within the context of the quasi-static approximation, have been studied in 
|52j . Using this approximation is very sensible because the hydrostatic time 
scale is very small, compared with stellar lifetimes, for many phases of the 
life of a star. It is of the order of 27 minutes for the sun, 4.5 seconds for a 
white dwarf and 10~^ seconds for a neutron star of one solar mass and 10 
km radius [53]. However, such an approximation does not apply to the very 
dynamic phases mentioned before. In those cases it is mandatory to take into 
account terms which describe departure from equilibrium, i.e. a full dynamic 
description has to be used. 

For the sake of generality, we have considered a locally anisotropic fluid. 
In fact, the assumption of local anisotropy of pressure, which seems to be 
very reasonable for describing the matter distribution under a variety of 
circumstances, has been proved to be very useful in the study of relativistic 
compact objects (see |l6l [51] and references therein). 

Finally we have also included viscous effects in our study. In fact, though 
they are generally excluded in general relativistic models of stars, they are 
known to play a very important role in the structure and evolution of neutron 
stars. Indeed, depending on the dominant process, the coefficient of shear 
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viscosity may be as large as rj ^ 10^° g cm~^ s~^ (see [55] for a review on 
shear viscosity in neutron stars). Also, a theorem by Raychaudhuri and De 
|56j . which states that in the evolution of non-dissipative charged dust the 
shear cannot vanish, emphazises the relevance of the shear in the evolution 
of charged fluids. 

On the other hand the coefficient of bulk viscosity may be as large as 10^° 
g cm~^ s~^ due to Urea processes in strange quark matter [57] (see also [5S] 
for a review on bulk viscosity in nuclear and quark matter). 

2 The energy- momentum tensor and the field 
equations 

In this section we provide a full description of the matter distribution, the line 
element, both inside and outside the fluid boundary, and the field equations 
this line element must satisfy. 

2.1 Interior spacetime 

We consider a spherically symmetric distribution of collapsing charged fiuid, 
bounded by a spherical surface S: we assume the fiuid to be locally anisotropic 
and undergoing dissipation in the form of heat fiow, free streaming radiation 
and shearing viscosity. For short we call this 'matter'. Choosing comoving 
coordinates inside S, the general interior metric can be written 

dsl = -A^de + B^dr^ + {Cr)\de^ + sin^ edcj)^), (1) 

where A, B and C are functions of t and r and are assumed positive. We 
number the coordinates = t, = r , x"^ = 6 and x^ = 0. 

The assumed matter energy momentum inside S has the form 

Taf3 = (/^ + P±)VaVf3 + P±gaf3 + (Pr " P±)XaXl3 + QaVp + V^g^ + eljp - 2r]aap, 

(2) 

where fi is the energy density, Pr the radial pressure, P± the tangential 
pressure, q" the heat fiux, e the radiation density, t] the coefficient of shear 
viscosity, V"' the four velocity of the fiuid, a unit four vector along the 
radial direction and a radial null four vector. These quantities satisfy 

i^"\4 = -i, \/"g„ = o, x"xa = i, x"K = o, rv; = -i, /% = o, (3) 
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and the shear a^p is given by 



(4) 



where the acceleration Oq, and the expansion O are given by 



(5) 



We do not exphcitly add bulk viscosity to the system because it can be 
absorbed into the radial and tangential pressures, and Pi, of the collapsing 
fluid p]. 

It should be noted that for a physically meaningful speciflc model we 
would need constitutive equations which would relate and determine the 
quantities /x, P^, Pi, g", e and rj. Without such relations we still have so 
many free functions that nothing useful can be said about the behaviour of an 
individual case. However, we will show that some important general physical 
results follow just from assuming, for example, that dissipation carries energy 
radially outwards. 

Since we assumed the metric ([1]) comoving then 



where g is a function of t and r. With ([6]) we obtain for (jl]) its non null 
components 



(Til = 



-P^CT, 0-22 



sin^^ 




(7) 



where 




(8) 



and the dot stands for differentiation with respect to t, which gives the scalar 
quantity 





(10) 



where the prime stands for r differentiation. 
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2.2 The electromagnetic energy tensor and the Maxwell 
equations 

The electromagnetic energy tensor is given by 

E-0 = ^ [f^'F^, - \F'''F,sga(^^ , (11) 

where Fa/s is the electromagnetic field tensor. Maxwell's equations can be 
written 

Fa/S = 4>l3,a - 4>a,p, (12) 

F'^f.p = Anr, (13) 

where 0q is the four potential and Jq is the four current. Since the charge is 
assumed to be at rest with respect to the coordinate system used in ([1]), there 
is no magnetic field present in this local coordinate system, and therefore we 
can write 

0„ = $(5°, J" = ?V^", (14) 

where the charge density, and $ are both functions of t and r. Charge 
conservation implies that 

s(r) =47r [\B{Crfdr, (15) 
Jo 

which is the electric charge interior to radius r, is time-independent. 
With ([1]) and ([6]) we obtain for the Maxwell equations (|T2l) and flT3l) 



B' 


-4 


'b 









'c 



(16) 
(17) 



Integrating f|T6|) and f|T7|) produces 



, sAB , , 

$ = . (18) 



7 



2.3 The Einstein equations 

Einstein's field equations for tlie interior spacetime ([1]) are given by 



al3 



The non null components of ^ with (fT|), (|2|), (|6|), (fTTj) and (fT8|) 

{sAY 



(19) 
become 



00) 



87r(Tol + Eoi) 
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22 



lly 



E. 



22) 



87r(/i + e)A^ + 



'B C\ C /A\ 



~r\B C) 
-8n{q + e)AB 
(C' _BC' 



C" 



' B' _ 

, 'b^'c 



'c 



21 



87r ( + e - ^r^a ) B' 



CA'\ 2 fB 
C^J ^r\B~C^ 
{sBf 



(20) 



(21) 



-V 

Aj 



'C 



c 



{CrY 
.AC 



- 2- 



AC 



(22) 



(C'Y A'C 









1 


(1 









in 



33 



E. 



33) 



^.(P. + ^w)(C.)^+(^)' 



B^C~A [b^c] ^BC 



A" C" A' fB' 
~A^~C~Ji\B 



'c 



B' C 

'b~c 



1 (A' B' 



The component (pT]) can be rewritten with ([8]) and (fTOl) as 
47r(g + 6)5 = ^(0- a)' -a^. 



(23) 



(24) 



8 



Next, the mass function m{t, r) introduced by Misner and Sharp (see 
also [61j) can be generahzed to include the electromagnetic contribution by 



m 



(Cr) 



D 23 



•-23 



Cr 
2 



rC 



{Cr)'' 
B 



n 2 



+ 1 



2Cr 



(25) 



which is the same mass function used in 



and [22]. 



2.4 The exterior spacetime and junction conditions 

Outside S we assume we have the Reissner-Nordstrom-Vaidya spacetime (i.e. 
we assume all outgoing radiation is massless), described by 



2M{v) 



dv^ - 2drdv + r^{de^ + r^ sin O^dcp'^) (26) 



where M(t>) and Q denote the total mass and charge respectively, and v is 
the retarded time. 

The junction conditions for the smooth matching of an adiabatic charged 
sphere to the Reissner-Nordstrom spacetime were discussed in detail in [22] 
and [^, whereas the matching of the full non-adiabatic sphere (including 
viscosity) to the Vaidya spacetime was discussed in [63]. The result is that 
the matching of ([T]) and fl26p on S implies 



Pr + 47]a = q, m{t,r) = M{v), s = Q, 



(27) 



where = means that both sides of the equation are evaluated on S. 



3 Dynamical equations 

The non trivial components of the Bianchi identities , (T~°^ + E~"f^)-p = 0, 
from (ED and 1^ yield 

{T--^ + E--^),pV^ = -l(^ + e)-(/. + P,, + 26-^r/a)-^ 

-2 (/i + P± + e + ^ ^ (28) 
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--(, + ey-2(, + e)l^ 

1 1 / 4 \' 

^W + ^) + ^(P. + ^-3w) 

+2(, + .)A + 2(, + e)iL 

/ A \ A' 



0, 



BCr AnB{Cr)* 



(29) 
0, 



where we used ( !T8|l and (fT5|) . 

To study the dynamical properties of the system, let us introduce, follow- 
ing Misner and Sharp [60], the proper time derivative Drp given by 



1 d 

Dt = 

^ Adf 



and the proper radial derivative Dr, 



where 



1 a 



R = Cr, 



(30) 



(31) 



(32) 



defines the proper radius of a spherical surface inside E, as measured from 
its area. 

Using (!30l) we can define the velocity U of the collapsing fluid as the 
variation of the proper radius with respect to proper time, i.e. 



U = rDj-C < (in the case of collapse). 
Then (1251) can be rewritten as 



E = 



(Cry 
B 



1 + U' 



2m{t, ''') _|_ ^ ■s 



Cr 



\CrJ 



1/2 



With ([3lD-([32]) we can express ^ as 

■1 



47r(g + e) = E 



a 



3/^.(6-.)-^ 



(33) 



(34) 



(35) 
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Using ([2nH2SD and ([3DH22D we obtain from 

4 



-Att 



Pr 



-7]a]U + {q + e)E 



and 



Dnin = Air 



/i + e + (g + e 



U 

'e 



(36) 



(37) 



Expression (136|1 describes the rate of variation of the total energy inside a 
surface of radius Cr. On the right hand side of (|36ll . (P^ + e — 4:r]a/3)U 
(in the case of collapse f/ < 0) increases the energy inside Cr through the 
rate of work being done by the "effective" radial pressure Pr — Arja/S and 
the radiation pressure e. Clearly here the heat flux q does not appear since 
there is no pressure associated with the diffusion process. The second term 
— (g + e)E is the matter energy leaving the spherical surface. 

Equation flTTI) shows how the total energy enclosed varies between neigh- 
boring spherical surfaces inside the fluid distribution. The first term on the 
right hand side of (137|) . /i + e, is due to the energy density of the fluid ele- 
ment plus the energy density of the null fluid describing dissipation in the 
free streaming approximation. The second term, {q + e)U/E is negative (in 
the case of collapse) and measures the outflow of heat and radiation. Finally 
the last term is the electrostatic contribution. 
Equation fl37|) may be integrated to obtain 

R r rr-1 o2 



m 



e+(g + e)- 



2R 2 



R s2 



(38) 



(assuming a regular centre to the distribution, so m(0) = 0). 

The acceleration Dj-U of an infalling particle inside S can be obtained 
by using ([22]), ([30]), ([25D and ([3lD, producing 

4 \ ^ EA' 
Pr + e r)a 



m 



DtU = - A-n 



i?2 



(39) 



and then, substituting A' / A from ([39]) into ([29]) . we obtain 



(/i + + 2e - fr/cr) DtU 
= - (/i + P^ + 2e - |r/a 



m 

+4:TT 



Pr + e- -r]a 



R ' 



-P2 



Dr (Pr + e- ha) +2 (Pr - Px + e 



-E 



2W)- 



47rP4 



Drs 



(40) 



Drq + Dre + 4(g + e)| + 2{q + e)a 
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which in the non-dissipative locally isotropic case coincides with Eq. (43) in 
[14J. Let us now analyze in some detail the three terms on the right of fl40p . 

The first term on the right hand side of PUI) represents the gravitational 
force. The factor within the round bracket (the same factor as on the left of 
fHOj) ) defines the inertial mass density ("passive" gravitational mass density) 
and shows how it is affected by dissipative terms. Observe that it is not 
affected by the electric charge. 

The factor within the first square bracket shows how dissipation and the 
electric charge affect the "active" gravitational mass term. Using (!38|) in (HO!) 
we see that the charge will increase the "active gravitational mass" only if 



otherwise it will decrease it. This strange effect was already noticed by 
Bekenstein [14], and enhances the possibility that Coulomb repulsion might 
prevent the gravitational collapse of the sphere. 

There are three different contributions in the second square bracket. The 
first one is just the gradient of the total "effective" radial pressure (which 
includes the radiation pressure and the influence of shear viscosity on P,,). 
The second contribution comes from the local anisotropy of pressure, includ- 
ing the contributions from the radiation pressure and shear viscosity. Finally 
the last term describes Coulomb repulsion, which is always positive (always 
opposing gravitation). 

The last square bracket contains different contributions due to dissipative 
processes. The third term within this bracket is positive {U < 0) showing that 
the outflow of g > and e > diminish the total energy inside the collapsing 
sphere, thereby reducing the rate of collapse. The last term describes an 
effect resulting from the coupling of the dissipative flux with the shear of the 
fluid. The effects of Dt€ have been discussed in detail in |6l]. Thus it only 
remains to analyse the effects of D^q] this depends on the transport equation 
adopted, and we will proceed to study one case in the next section. 

However before doing that it is instructive to recover a known result for 
the static case. 




(41) 



or, equivalently 




(42) 
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3.1 Static charged dust 



In the limit of hydrostatic equihbrium when U = a = q = e = 0, we have 
from ([29]) 

P; + (, + p4 + 2(P.-Pj(^-j^^0. (43) 

which is just the generahzation of the Tolman-Oppenheimer-Volkov equa- 
tion for anisotropic charged fluids obtained in [62j while studying dynamical 
instability for radiating anisotropic collapse. 

When the static fluid reduces to charged dust, with Pr = P± = 0, then, 
by using f[T51) . fH5]) becomes 

Since B and C depend only on r, we can transform r so that B = C. 
Eliminating s from the field equations fl22p and (1231) . we can solve for AB, 
and imposing regularity conditions and rescaling t we have 

C = B, AB = 1, = r^B'\ (45) 

Substituting into flUl) we obtain 

fx' = (46) 

which is the well known result originally obtained by Bonnor [5] for arbitrary 
symmetry. 



4 The transport equation 

We shall use a transport equation derived from the Miiller-Israel-Stewart 
second order phenomenological theory for dissipative fluids [651 [66] . 

Indeed, it is well known that the Maxwell-Fourier law for radiation flux 
leads to a parabohc equation (diffusion equation) which predicts propagation 
of perturbations with infinite speed (see [67]- [70] and references therein). 
This simple fact is at the origin of the pathologies [7T] found in the ap- 
proaches of Eckart [72] and Landau [73] for relativistic dissipative processes. 
To overcome such difficulties, various relativistic theories with non-vanishing 
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relaxation times have been proposed in the past ESI EH [75]. The impor- 
tant point is that all these theories provide a heat transport equation which 
is not of Maxwell- Fourier type but of Cattaneo type [76], leading thereby to 
a hyperbolic equation for the propagation of thermal perturbations. 
The corresponding transport equation for the heat flux reads 



(47) 



where h'^'^ is the projector onto the three space orthogonal to V^, k denotes 
the thermal conductivity, and T and r denote temperature and relaxation 
time respectively. Observe that, due to the symmetry of the problem, equa- 
tion fH71) only has one independent component, which may be written after 
using m. P and (HM as 



rq 



K 

B 



(TA)' - qA. 



(48) 



Now using (I30tl34p we can rewrite ( l48l) as 



i^T^qr^ f r \ /3U 1 1 



kE ^ kT 
—DrT - —DtU 

T tE 



1 



(49) 



We can couple the transport equation in the form above, ( l49l) . to the 
dynamical equation (HOj) . in order to bring out the effects of dissipation on 
the dynamics of the collapsing sphere. For that purpose, let us substitute 
( 149|) into (HOj) : then we obtain, after some rearrangements. 



/i + + 2e - -ar]j (1 - a)DTU 

kE'^ 

= {l-a)Fgrav + Fhyd+ DrT + E 



-Eq 



5U 
2B 



-a 



2Ee 2 



KT^q 
2r 

U 
> 

'B 



Dn 



kT^ 



-Dre 



a 



(50) 



where F„rav and Fhyd are defined by 



14 



grav 




-E^ \^Dr [Pr + t- -r^Or) 

+2(P, -Px + e-2r/or)4- 



R P2 



(51) 



P, 



hyd 



S 



DrS 



(52) 



47rP4 



and a is given by 




— U + P^ + 2e - -ar] 



-1 



a = 



(53) 



Some comments are in order at this point: 

• Once the transport equation has been taken into account, then the in- 
ertial energy density and the "passive gravitational mass density", i.e 
the factor muhiplying DtU and the first factor at the right of fHUl) 
respectively (which of course are the same, as expected from the equiv- 
alence principle), appear diminished by the factor 1 — a, a result already 
obtained in [15], but here generalized by the inclusion of the viscosity 
and radiative phenomena. 

• Observe that the charge does not enter into the definition of a. However 
it does affect the "active gravitational mass" (the factor within the 
square bracket in flSTl) ). 

• The repulsive Coulomb term (the last term in 0521] ) depends on Dpts 
and always opposes gravitation. Its effect is reinforced if Drs is large 
enough to violate (H2|) . in which case the charge will decrease the "active 
gravitational mass" term in (ISTj) . 



In this section we shall find some interesting relationships linking the Weyl 
tensor with matter variables, from which we shall extract some conclusions 
about the arrow of time. 
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The Weyl tensor 
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From the Weyl tensor we may construct the Weyl scalar = C^^^C^p^s 
which can be given in terms of the Kretchman scalar IZ = R^^'^^ Rap'^s, the 
Ricci tensor Rap and the curvature scalar R by 



(54) 



Substituting (172|) from the Appendix with the field equations ( !20ll23|) into 
(!54l) we obtain 



477- 

8 = m- —{fx -P, + P^ + 2r]a)R' - -, 



where £ is given by 



S 



C 



481/2 



R' 



(55) 



(56) 



From fl55l) with (l36l) and fl37j) we have 
■1 



fi + P^ + e+^r]a) R^U + {q + t)ER^ 



+ ^. (57) 



and 



Dr£ = 4tt 



(g + e)^ - ^R'Dnifi -Pr + P± + 2r]a) 



+ (e + Pr -P±- 2ria)R^ 



R \R 



(58) 



From (1581) we obtain at once for the non-charged, non-dissipative, perfect 
fluid case 

Dr£ + —R'Dnfx = 0, (59) 

implying that Dr^i = produces C = (using the regular axis condition), 
and conversely the conformally fiat condition implies homogeneity in the 
energy density. 

This particularly simple relation between the Weyl tensor and density 
inhomogeneity, for perfect fluids, is at the origin of Penrose's proposal to 
provide a gravitational arrow of time in terms of the Weyl tensor [77]. The 
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rationale behind this idea is that tidal forces tend to make the gravitating 
fluid more inhomogeneous as the evolution proceeds, thereby indicating the 
sense of time. 

However the fact that such a relationship is no longer valid in the pres- 
ence of local anisotropy of the pressure and/or dissipative processes, already 
discussed in [IB] , explains its failure in scenarios where the above-mentioned 
factors are present. Here we see how the electric charge distribution affects 
the link between the Weyl tensor and density inhomogeneity, suggesting that 
electric charge (whenever present) should enter into any deflnition of a grav- 
itational arrow of time. 

6 Conclusions 

We have provided a full set of the equations required for a description of phys- 
ically meaningful models of collapsing charged spheres. We have included 
dissipative phenomena as well as anisotropic pressure; the justiflcation for 
doing so was given in the Introduction. 

The role of charge distribution in the dynamics of such conflgurations 
is clearly exhibited in equations fl36|) . fl37|) . fHOl) and fl50|) . In particular it 
is worth stressing the fact that electric charge, unlike pressure, does not 
always produce a "regeneration effect" (does not always increase the "active 
gravitational mass"). This fact together with the presence of the Coulomb 
term in (l40l) (or (l50l) ) indicates the relevance of the electric charge in the 
process of collapse. 

Finally we have obtained a relation (158!) exhibiting the way in which 
electric charge affects the link between the Weyl tensor and density inhomo- 
geneity. The consequences of this for a definition of a gravitational arrow of 
time have been discussed. 
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Appendix 



The spacetime ([T]) has the following non null Riemann tensor components 



-Roioi 
-R0202 

-R0212 

-R1212 

-R2323 



-BB + -Ab + AA" - -A'B', 
A B 



{Crf 
{Crf 
{Crf 

{Crf sin^ 6 



C AC /AVA'fC V 

c^Ac^ \b) 'Jy'c^r, 



C BC C A! 1 (B C' 
'C ^ B~C ^ 'C~A ^~r\B^~C . 



a) BC 

'rC\ 
C I 



C" 


B' C 


~c 


^ 'b~c 






-( 











-4)1 



CrC 



C\^ , 
b) +^ 



(60) 
(61) 

(62) 

(63) 

(64) 



and 

-R0303 = -R0202 sin^ -R0313 = -R0212 sin^ 6', -R1313 = -R1212 sin^ (65) 

hence it has 5 independent components and the Kretchman scalar becomes 

r 1 2 4 

^ = ^ [(AB)^^^°^°^)' + (IC;0^^^°^°^)' ' \AmCr)'^^"""' 

2 2 1 2 

"'"TrTFm (-^1212) + NO . 4 ^ (-R2323) 

The components ( l60ll64l) can be written in terms of the Einstein 
— 5'a/3R-/2 and the mass function ( l25l) producing 



(66) 
tensor G^/j = 



-Roioi = {ABf 
-R0202 = 
-R0212 = 
-R1212 = 
-R2323 = 



2^2^00 - 21^^11 + (^^22 - (m - 



(5Cr)2 



1 

(Cr)3 



2Cr) 



2^ ^00 ~ 

2Crsin2^(m- ^ 



[m — 

2C> 



2Cr, 



(67) 
(68) 
(69) 
(70) 
(71) 
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Substituting ([SMID into 
48 / 



n 



we obtain 
16 / 



(Cr) 



m 



2Cr 



(Cr) 



m 



2Cr 



00 



\ab) 



+ 3 



GooGn 



00 



G 



11 



52 
G22 



Goo Gn 
'A^ ~ 'B^J (Cr)2' 



+ 4 



G 



A^ 52 
2 



Cll _|_ ^22 



(Cr) 



22 



(Cr) 



(72) 
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